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In a network of nonlocally coupled Stuart-Landau oscillators with symmetry-breaking coupling, 
we study numerically, and explain analytically, a family of inhomogeneous steady states (oscillation 
death). They exhibit multi-cluster patterns, depending on the cluster distribution prescribed by 
the initial conditions. Besides stable oscillation death, we also find a regime of long transients 
asymptotically approaching synchronized oscillations. To explain these phenomena analytically in 
dependence on the coupling range and the coupling strength, we first use a mean-field approximation 
which works well for large coupling ranges but fails for coupling ranges which are small compared to 
the cluster size. Going beyond standard mean-field theory, we predict the boundaries of the different 
stability regimes as well as the transient times analytically in excellent agreement with numerical 
results. 
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I. INTRODUCTION 

Coupled nonlinear systems exhibit a plethora of differ¬ 
ent collective behavior such as partial or cluster synchro¬ 
nization dH, spatio-temporal patterns under delayed 
feedback control dH, and chimera states which consist 
of spatially coexisting domains of synchronized and un¬ 
synchronized dynamics [TOHE]- Recently, special atten¬ 
tion has been paid to the different types of oscillation 
quenching, amplitude death and oscillation death nsr 
ng. Amplitude death and oscillation death differ in the 
mechanisms by which they are induced: In coupled os¬ 
cillator networks, the coupling between oscillators can 
stabilize an already existing homogeneous steady state 
which is unstable in the absence of coupling. This phe¬ 
nomenon is called amplitude death. In contrast, oseilla- 
tion death requires the breaking of the system’s symme¬ 
try: Here inhomogeneous steady states are newly created 
as well as stabilized through a symmetry-breaking cou¬ 
pling between the oscillators. A combination of these two 
phenomena of spontaneous symmetry-breaking, chimera 
states and oscillation death, has recently been found in 
a paradigmatic model of nonlocally coupled supercritical 
Hopf normal forms (Stuart-Landau oscillators) pTl [TS] . 
The resulting patterns of coexisting coherent and inco¬ 
herent inhomogeneous steady states occur in form of clus¬ 
ters of different size, and have been named ehimera death. 
Similar chimera death patterns have recently also been 
observed with global (all-to-all) coupling [T9j [20] . 

Oscillation death has been observed experimentally in 
many different systems, such as chemical reactors [21], 
chemical oscillators [22] , chemical droplets [23] , electronic 
circuits [24] , or thermokinetic oscillators [25] . There exist 
various biological applications, e.g. neural networks [26] . 
genetic oscillators [27], calcium oscillators [28] or stem 
cell differentiation [29] where it has been proposed as a 
basic mechanism for morphogenesis and cellular differ¬ 
entiation. Consequently, there has also been theoretical 
interest in oscillation death and it has been shown to 


exist for special time-delayed m and repulsive [3I1I32] 
types of coupling as well as coupling through conjugate 
or dissimilar variables [33ll35] . However, research has 
mostly focussed on small numbers of coupled oscillators 
in networks with local or global coupling. 

The influence of nonlocal coupling on the onset of oscil¬ 
lation death and on the emergence of multi-cluster pat¬ 
terns in larger networks has not been investigated sys¬ 
tematically, so far. It is possible to use the coupling range 
P as a control parameter, e.g., in a bidirectionally cou¬ 
pled ring of N nodes. Thus one can interpolate between 
the well-studied limit cases of local (nearest neighbor, 
P = 1) and global (all-to-all, P = [A^/2j) coupling. 

In our present study, we analyze the transition 
from completely synchronized oscillations to multi¬ 
cluster inhomogeneous steady state patterns (oscillation 
death) in networks of nonlocally coupled oscillators with 
symmetry-breaking coupling. We choose initial condi¬ 
tions which favor the emergence of spatially coherent 
clusters, in contrast to m, where hybrid patterns with 
coexisting domains of spatially coherent and spatially in¬ 
coherent steady states (chimera death) were analyzed. 
We investigate the effects of several factors on the oc¬ 
currence of different oscillation death states both numer¬ 
ically and analytically. First, we use nonlocal coupling 
and thus cover the whole range between the limit cases 
of local and global coupling. Second, we also vary the 
initial conditions, specified by cluster distributions, to 
find different coexisting spatial patterns for oscillation 
death. Furthermore, we also vary the coupling strength, 
and combine all three variations. Beside the regions of 
stability, we also find parameter regimes where the os¬ 
cillation death state is transient, but persists for a long 
time. Due to the coherent nature of the oscillation death 
patterns, we are able to predict analytically and with 
high precision the boundaries between different stability 
regimes. 

This paper is organized as follows: In section |H] we in¬ 
troduce our model of nonlocally coupled oscillators and 
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briefly summarize previous work on two coupled oscilla¬ 
tors m- Our approach is numerical in the first step, and 
we observe two different types of oscillation death: tran¬ 


sient as well as asymptotically stable. In section III 


present our numerical findings. In the next step, we ex¬ 
plain our findings analytically, such as boundaries in the 
parameter space for the occurrence of the different pat¬ 
terns of oscillation death. In section [IV| we use a mean- 
field ansatz to obtain a first approximation of the stabil¬ 
ity boundaries. In sect ion |V| we go beyond the mean-field 
approximation to describe the boundaries analytically in 
excellent agreement with the numerical results. In sec¬ 
tion [VT] we calculate the transient times of transient oscil¬ 
lation death analytically, thereby gaining insight into the 
mechanism. In section EIIl we draw some conclusions. 


II. MODEL 

In the present study, we investigate a network of cou¬ 
pled Stuart-Landau oscillators. The equation of the sin¬ 
gle Stuart-Landau oscillator, which is a truncated generic 
expansion of a system near supercritical Hopf bifurcation 
in center manifold coordinates, is given by 

Z = f{z) = {X + iLO - \z\^) z. (1) 

Here z = = x^iy e C denotes the phase space, and 

A G M can be viewed as a bifurcation parameter. A Hopf 
bifurcation occurs at A = 0. For A > 0, a stable periodic 
orbit with radius vq = \/A appears, which oscillates with 
angular velocity cj G M. This periodic orbit z{t) = 
is harmonic and rotationally symmetric. In mathemat¬ 
ical terms, we say that the single Stuart-Landau oscil¬ 
lator is S'^-equivariant, i.e. z{t) solves i = f{z) if and 
only if e^^z{t) solves i = /(^), for any fixed phase angle 
0 G [0, 27r] = S^. Indeed / commutes with the S^-action 
of 0 on C^: f{e^^z) = f{z) for all 0 and 

We arrange N Stuart-Landau oscillators Zj^ j = 
1, 2,..., A", with N even, as a ring and couple them as 
follows: 


= fi^j) + ^ Yj Re Zj) (2) 

k=j-P 

where all indices are modulo N. The coupling strength is 
described by the real parameter a > 0. Each oscillator 
Zj is coupled to the oscillators Zj-i ,..., Zj-p as well as 
to the oscillators ..., i.e., to its P nearest 

neighbors in each direction of the ring, respectively. The 
coupling range is given by P, where P corresponds to the 
number of coupling neighbors in each direction on a ring. 
The coupling is normalized by the number of existing 
links, which is 2P. 

In this publication, we focus on the case of nonlocal 
coupling. Note, however, that our results also hold true 
for the case of local coupling, which is given by P = 1, 
and in the limit of large N for global coupling, which 


is given approximately by P = A/2. Strictly speaking, 
the limit of global coupling exists only for odd A and 
P = (A —1)/2, since for even A the antipodal link would 
appear twice. 

Here we consider coupling only in the real part Xj of 
the variable zj^ which breaks the P^-equivariance of the 
single Stuart-Landau oscillator: Re{P'^Zj) = Re{zj) 
is only true for ^ = 0 or ^ = tt. Symmetry-breaking of 
P^-equivariance is a necessary condition for the existence 
of isolated nontrivial steady states with Zj ^ 0 and thus 
for oscillation death. 

Note, however, that the rotationally symmetric peri¬ 
odic orbit of the single Stuart-Landau oscillator is pre¬ 
served by the coupling. It exists for all coupling strengths 
(j and can be seen if all oscillators are in synchrony. 

In the present study, we restrict our parameters to the 
interesting case 0 < A < cj. The following results, which 
will be needed later, have been derived in [30] for a system 
of only two coupled Stuart-Landau oscillators: 

Zi= f{zi)+e{Rez 2 -Rezi) 

Z 2 = fiz 2 ) + £ (Re ^1 - Re Z 2 ) , 

with a coupling parameter e > 0. In the in-phase 
subspace zi = 2 ( 2 , this system simplifies to i = /(^), 
i.e., to a single Stuart-Landau oscillator. In the anti¬ 
phase subspace zi = — 2 ( 2 , this system simplifies to 
i = f{z) — 2sliez. The system i) exhibits a stable 
synchronous periodic orbit for all coupling parameters 
5 . There also exists a trivial homogeneous steady state 

= 2(2 =0. Linearization at zero reveals that this state 
is always unstable. Supercritical pitchfork bifurcation 
occurs at 

ep = 5 (a + ^) . (4) 

At this pitchfork bifurcation, two inhomogeneous steady 
state solutions {xi^ 2 ,yi, 2 ) 

Xi ^2 = ±\/ X-e-\-\/e‘^- up-y‘^ (5) 

2/1.2 = (As - _|_ Av^ - / ( P ) ( 6 ) 

emanate for e > ep. They correspond to antiphase states 
zi = —Z 2 . The radius f of the inhomogeneous steady 
state is given explicitly by 

f^ = A — £ -h \/ (7) 

Linearizing at those steady states, we find the eigenvalues 
T]. There exists a pair of complex conjugate eigenvalues 
in direction of the in-phase subspace, 

Ml ,2 = A - 2f^ ± i\/cc;2 (8) 

and there exist two real and distinct eigenvalues in direc¬ 
tion of the anti-phase subspace, 

Ms, 4 = —A + £ — 2^/± (A — 5 ). (9) 
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Figure 1. Top: Real part of the complex conjugate eigenval¬ 
ues (green dashed) and of the real eigenvalues and 
(gray solid) as a function of the coupling strength £ for 
2 coupled oscillators. Bottom: Imaginary variable y of the 
inhomogeneous steady states versus s. Parameters: A = 1, 
uj — 2. There is a pitchfork bifurcation (PB) oX e — 2.5, in¬ 
dicated by a diamond, and secondary Hopf bifurcations (HB) 
at £ = 4.25, indicated by circles. 


In particular we encounter secondary Hopf bifurcations 
at 

effB = j(A + 4^), (10) 

where the inhomogeneous steady states are stabilized, 
thus marking the onset of oscillation death. The sec¬ 
ondary eigenvalues /ii ^2 remain complex conjugate eigen¬ 
values, with relatively small real part. In contrast, /i 3^4 
give two distinct real and negative eigenvalues. See Fig.j^ 
for a visualization of the eigenvalues and inhomogeneous 
steady states. 

III. NUMERICAL SIMULATIONS 

As a first step in the analysis of the system, we solve 
Eqs. numerically. We vary the coupling range P 
from 0 (no coupling) to N/2 (global coupling) in steps 
of 1. The intermediate range 0 < P < N/2 is called 
nonlocal coupling and our main interest lies here. Ad¬ 
ditionally, we also increase the coupling strength a from 
cr = 1 to 30 with stepsize 1. The simulations are run 
no longer than time t = 4000. The other parameters are 
fixed to N = 100, A = 1, and uo = 2. As initial con¬ 
ditions, we use anti-phase clusters of size n: In a first 
step, all oscillators zi,... ,Zn of the first cluster are set 
to the values (xj^yj) = (—1,+1), which we call the up¬ 
per branch since yj is close to the upper steady state y 


(cf. bottom panel of Fig. [^. In contrast, the oscilla¬ 
tors 2 )^+ 1 ,..., 2 : 2 n are set to values on the lower branch 
with {xj^yj) = —( —1,+1) = (+1,-1). Then oscilla¬ 
tors Z2n-\-i^ • • • 5 ^3n are again on the upper branch and 
so forth. In a second step, a random number, drawn 
form an underlying Gaussian distribution with zero mean 
and variance 0.1, is added to the value of each os¬ 
cillator. We have also performed simulations with dif¬ 
ferent initial conditions, such as (xj^yj) = (0.2,—1) or 
{xj^yj) = (—v^O^, Ty/^). We have not found any no¬ 
ticeable difference in the asymptotic behavior for those 
initial conditions. The simulations are performed for 
cluster sizes n = 50, n = 25, n = 10, n = 5, n = 2, 
and n = 1. 

The asymptotic results are depicted in Fig. in the 
(P,cr)-plane for different initial conditions (panels a - e), 
exhibiting high multistability. As discussed above, the 
initial conditions are regular anti-phase clusters of dif¬ 
ferent size with superimposed random additions. For 
those initial conditions amplitude chimera states cannot 
be observed, in contrast to mi m where other initial 
conditions were used and also transients were recorded. 
We observe a large region of synchronous oscillations for 
each coupling range P. We also find regions of sta¬ 
ble oscillation death, color-coded by yellow (light), or¬ 
ange (medium), and red (dark), depending on the cluster 
size. By m-cluster oscillation death (m-OD), we denote 
a steady state such that there exist m clusters on the 
upper branch as well as m clusters on the lower branch. 
These clusters do not necessarily have the same size. We 
observe a large region of oscillation death where the clus¬ 
tering takes the same form as given by the initial condi¬ 
tions. This region can in general be found for large cou¬ 
pling strength a and large coupling range P/N. It is the 
boundary of this specific region that we investigate ana¬ 
lytically below. Furthermore, we also find regions where 
there are asymptotically three times as many clusters as 
in the initial conditions (e.g., 3-OD in panel (a), 6-OD 
in panel (b), 15-OD in panel (c) of Fig. [^, and regions 
where we find five times as many (e.g., 5-OD in panel 
(a), 10-OD in panel (b)). These regions are marked by 
darker color. A scenario leading from 1-OD to 3-OD is 
shown in Fig. In panel (b) the 1-OD state is unstable, 
and the two oscillators in the middle of the cluster switch 
to the opposite branch. In panel (c) two more oscillators 
switch to the opposite branch. 

We observe that, within one cluster, not all oscillators 
assume the same values for radius and phase: We call 
this effect cluster deformation (Fig. |^. The deforma¬ 
tion is particularly pronounced for small coupling range 
(Fig. [^), and is approximately linear, either decreasing 
or increasing, seen from the corner oscillators of the clus¬ 
ter. In the middle of the cluster, we observe a plateau, 
where the oscillators have the same radius and phase. We 
will describe this phenomenon in more detail in section 
|V| where it will be used to find a better approximation 
of the behavior of system (§ beyond mean-field theory. 

Additionally, in the region where the asymptotic be- 
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Figure 3. Transition from 1-OD to 3-OD. (a) P = 26. (b) 
P = 25. (c) P = 24. Parameters: N = 100, A = 1, Ct; = 
2, cr = 25, n = 50. The simulations were run up to time 
t = 500. 



Figure 4. Cluster deformation effects of oscillation death for a 
coupling strength cr = 15 and different values of the coupling 
range P/N. Snapshots for the variables yj. (a) P — 10, (b) 
P = 20, (c) P = 25. Other parameters: N = 100, A = 1, 
u — 2^ n — 25, m — 2. 


Figure 2. Numerical simulation of the asymptotic behavior 
of the system in dependence on the coupling range P 
(horizontal axisj^and the coupling strength a (vertical axis). 
The green (gray) region marks synchronized oscillations, the 
colors yellow (light) to red (dark) mark the various m-cluster 
oscillation death (m-OD) states. White patches correspond to 
more complex patterns. The insets show the asymptotic OD 
state corresponding to P = 35 and cr = 20 (marked by white 
squares). Initial conditions: (a) n — 50, (b) n — 25, (c) n — 
10, (d) n = 5, (e) n = 2 (cluster size n = N/{2m)). Analytical 
approximations: White diamonds: mean-field approximation 
ap^n- Light blue circles: beyond mean-field approximation 
(Tpn- Parameters: N = 100, A = 1, cj = 2. 


havior is given by a synchronized oscillation (SYNC), we 
find transient oscillation death. Fig. depicts the time 
series of j = 1,..., 100, in the transient regime. It 
shows a transition from the initial state to transient 3- 
OD. For a long time, the inhomogeneous steady states 
seem stable but, eventually, a second transition to syn¬ 
chronous oscillations occurs. Following these observa¬ 


tions, we have also recorded the transient times T, see 
Fig-E We generally observe an increase of T if the cou¬ 
pling strength a is increased. The transient time will be 


investigated in more detail in section VI 


IV. MEAN-FIELD APPROXIMATION 

It is the purpose of this section to introduce a simpli¬ 
fied model, which is able to explain the above numerical 
observations. In this simplified model, we assume that 
we have clusters of size n, where the single components 
of the cluster are supercritical Stuart-Landau oscillators 
Eq. 0 - Since we restrict the considerations to cluster 
distributions with uniform cluster size n, we obtain the 
condition N = 0 mod 2n. For N = 100, e.g., this is the 
case for n = 50, 25,10, 5, 2 and I. We assume P-nearest 
neighbor coupling, see Eq. (2). 

As a simplifying mean-field assumption, and deviating 
from the numerics, we assume in this section, and only in 
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Figure 5. Time series yj of transient oscillation death: the 
hrst transition is from the initial condition to an unstable in¬ 
homogeneous stationary state that shows a second transition 
to synchronous oscillations, i.e., transient oscillation death. 
Parameters: N — 100, A = 1, Ct; = 2, cr = 6, P = 30. 
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Figure 6. Transient time T in the plane of the coupling range 
P/N and the coupling strength a with 1-OD (n = 50) initial 
distribution. The colorcode (grayscale) shows the transient 
time T ; the light gray region denotes stable oscillation death. 
Parameters: A/’=100,A=l,Ct; = 2. As initial conditions we 
used the explicit values given by Eqs. and (|^, but substi¬ 
tuted the coupling strength a instead of e. Here, no additional 
random numbers were added to the initial conditions. 
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Coupling range P 


this section, that all oscillators on the upper branch have 
the same radius and phase. We also assume that they 
are anti-symmetric (anti-phase) with respect to the os¬ 
cillators on the lower branch, implying that those also all 
have the same radius and phase. Therefore the coupling 
term vanishes within the clusters. It is between the clus¬ 
ters that the coupling term gives non-zero contributions 
(coupling of the real parts, inducing the inhomogeneous 
steady states, cf. section 0 - 

Our goal is to replace the model of N coupled oscil¬ 
lators by a simplified model of two Stuart-Landau os- 



Figure 7. Scheme of nonlocal coupling. For the oscillator j 
and P = 4 we have 4 relevant links to the left (solid red) as 
well as 4 irrelevant links to the right (dashed green). Irrele¬ 
vant links do not contribute to the coupling term. 


dilators which are coupled via their real parts with an 
effective coupling strength 


Zi = f{zi) + Kp^nCr {Rez 2 -Rezi) 

-22 = f{z 2 ) + Kp^nO- (Re Zi - Re 2 : 2 ) 

Since the dynamics of this system can be described an¬ 
alytically, cf. section [nj we can draw conclusions about 
the onset of oscillation death, provided that we know the 
effective coupling strength Our task hence is to 

find the scaling factor np^ni depending on both the cou¬ 
pling range P and the cluster size n, as indicated by the 
numerical simulations. 

As a first reduction step, note that it is sufficient to 
consider only the links emanating from one cluster in¬ 
stead of those of the complete network. Indeed the net¬ 
work structure implies that the coupling structure is the 
same for every cluster. This reduces the number of oscil¬ 
lators for which we have to calculate the effective coupling 
from N to n. 

In the following, we number the oscillators within one 
cluster with indices j = 0 ,...,n — 1, where n is the size of 
the cluster. The oscillator with index j is coupled to 2P 
neighbors: P to the left, and P to the right. We count 
^relevant’ links (+1), which are those between oscillators 
of different clusters. In contrast, irrelevant’ links, which 
are the links between oscillators of the same cluster, are 
not counted (0); see Fig. 

As a first step, consider a coupling range P < 2n. For 
the oscillator j we now count P links to the left. Of these, 
a maximum oi P — j is relevant, since the first j links 
to the left are between oscillators of the same cluster. 
If P — j is negative, then all left links are irrelevant. 
The number of relevant left links is therefore bounded 
from below by zero. The upper bound of the number 
of relevant links to the left is n, corresponding to the 
full cluster being relevant. In conclusion, the number of 
relevant links to the left of oscillator number j is given 

by 


max { min {n, P — j} , O}. (12) 

Now also consider coupling ranges P > 2n, more pre¬ 
cisely, consider P G {2nk + 1, • • •, 2n{k + 1)}, k G Nq. 
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Here k = 0 corresponds to P < 2n. Note that in this 
case the respective lower bound of relevant links to the 
left is given by kn ( i.e. k full clusters), while the upper 
bound is given by {k-\-l)n (i.e. {k-\-l) full clusters). Con¬ 
sequently, for the oscillator j, we find that the number of 
relevant links to the left is given by 

max { min {(/c + 1) n, P — j — kn} , kn}. (13) 

By reflection, the number of relevant links to the left of 
oscillator j is the same as the number of relevant links to 
the right of oscillator n — j — 1. 

In the next step, in order to reduce from n oscillators to 
two oscillators, we calculate the mean relevant eoupling 
of the whole cluster, and therefore consider each cluster 
as lumped into only one effective oscillator. We sum up 
the respective value of each oscillator within the cluster, 
and divide by the number of neighbors P as well as by the 
number of oscillators n to obtain a first approximation 
for the average coupling npy. 




Pn 


n—l 

max { 

j=0 


min {(/c + 1) n, P — j — kn} , kn} 


(14) 

This can be inserted into the threshold condition for sta¬ 
ble oscillation death obtained in section [h} shb = i(^ + 
4cj^/A). In total, for P G {2A:n + 1,..., 2 (/c + 1) n}, 
/c G No, we find the threshold 


_ _ I (A + 4u;VA) Pn _ 

{ min {(/c -h l)n, P — j — kn} , kn} 

(15) 

for the onset of oscillation death. 

Note that the size of the system, i.e., the total number 
of oscillators does not appear in the threshold. 

We can compare this analytic mean-field result with 
the numerical simulations: For a coupling range P which 
is of the same order of magnitude as the cluster size n, 
or larger, the analytic and the numeric thresholds agree 
very well. For a coupling range P which is small com¬ 
pared to the cluster size, we note that the analytic results 
give smaller thresholds than the simulations, see white 
diamonds in Fig. This suggests that the simplified 
mean-field model must be improved for this parameter 
regime. 


V. BEYOND MEAN-FIELD THEORY 

As mentioned above, the assumption that all oscilla¬ 
tors within a cluster have the same radius and phase is 
not true in general and deviates from the numerical re¬ 
sults. Therefore, we aim for a correction term to the 
mean-field theory used above, depending on the coupling 
range P and the cluster size n. 

In our reduced model, we replace the mean- 
field coupling by the corrected mean coupling 


{i^p,n — i^pn) section to es¬ 

timate the correction term Bp^. 

In Fig.|^we observed a deformation of the cluster shape 
(branch splitting) for small P. The deformation is ap¬ 
proximately linear, where the x- and the ^-variable of the 
two corner oscillators of the cluster stay (approximately) 
the same. The values of the neighboring oscillators, pro¬ 
ceeding inwards, either decrease or increase. Seen from 
the left, there are min{P mod n, (P 1) mod n, n — (P + 
I)modn,n — Pmodn} oscillators which show a linear 
decrease/increase with respect to their left neighbor. To 
first order, it follows that there exists a plateau, where 
the oscillators have the same radius and phase. On the 
right we encounter an analogous linear increase or de¬ 
crease, respectively. The cluster is reflection symmetric 
in the sense that the phenomenon is the same if we ap¬ 
proach from the left or the right. 

This behavior is due to the fact that not all oscillators 
“feel” the same coupling: The corner oscillators feel more 
or less oscillators from the other clusters than those in 
the middle, depending on the precise values of P and n. 

An exception are the values P G {nk, nk — 1 : k e N} 
where the mean coupling coincides exactly with the cou¬ 
pling that every single oscillator feels; thus no branch 
splitting effects occur (Fig. |^). In particular, we can 
conclude that, for those P, no correction term is neces¬ 
sary. 

Following the description above, the effect of the clus¬ 
ter deformation on the mean coupling is proportional 
to min{P mod n, (P + I) mod n, n — (P -1- I) mod n, n — 
Pmodn} as well as proportional to the cluster size n. 
Overall, the significance of the correction term decreases 
as P increases, since the more oscillators are coupled, the 
closer we are to the mean-field approximation. A simple 
intuitive ansatz is to assume a global 1/P dependence of 
the correction term. Thus, to compensate the cluster de¬ 
formation term, we need a prefactor P~^ to achieve the 
desired global I/P dependence. 

We suggest the following term: 

Bpn — cnP~^ min{Pmodn, (P -j- I) modn, 

n — (P -j- I) mod n^n — P modn}. 

( 16 ) 

Here, we denote the proportionality constant by c. We 
have found it to be the same for all our examples, 
c ~ 0.047. We therefore assume that this constant is in¬ 
dependent of n and P. Furthermore, we also claim that 
it is independent of A, since the single oscillators are not 
influenced by the circumference of the ring. However, 
there might be an additional dependence on A and cj, 
which is not the subject of our present investigation. 

As a consequence, we also obtain a corrected formula 
for the oscillation death threshold, which we denote by 


^P,n 


A + 4ci;^/A 
4(/^P,n ~ ^P,n) 


(17) 
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Note that Eq. 0 gives either a larger threshold or the 
same threshold as before. The difference of the simpli¬ 
fied and the improved formula is most noticeable if P 
is small compared to the cluster size n, compare Fig. 
The improved formula gives excellent agreement with the 
numerical threshold for stable oscillation death. 


VI. TRANSIENT BEHAVIOR 

The aim of this section is to describe the tran¬ 
sient times analytically. Transient oscillation death, see 
e.g. the time series depicted in Fig.[^ occurs for coupling 
strengths a between the emergence of the unstable inho¬ 
mogeneous steady state branches (via a pitchfork bifur¬ 
cation, cf. section [n| and [30]) and the stabilizing Hopf bi¬ 
furcation. The manifestation of the long transient times 
is due to different orders of magnitude of the real parts 
of the eigenvalues 

Re/ii^2 = -A + 2Rp^nO' - (18) 

Re/Us = 1^3 = (19) 



Figure 8. Transient time T vs. coupling strength a for dif¬ 
ferent variances of the initial distribution for 1-OD. Initial 
conditions {xj,yj) = (—0.2,-h0.7) on the upper branch, and 
(xj^i/j) = (-1-0.2,—0.7) on the lower branch, with random 
fluctuations of different variance added. Green (lighter) dots: 
Numerical values using a variance of 10~^. The das hed green 
line gives the analytical values calculated in Eq. ( [^ with 
fit parameter {3 — 1.22. Magenta (darker) dots: Numerical 
values using a variance of 10“^. T he s olid red line gives the 
analytical values calculated in Eq. (j^ with (3 — 0.382. Pa¬ 
rameters: N — 100, A = 1, ct; = 2, Rp^n — 0.48898, P — 48, 
n = 50, m = 1. 


Here we have used the abbreviation Rp^n = i^p,n ~ ^Pn* 
The absolute value of jxs increases with increasing cou¬ 
pling, see also Fig. for the system of two coupled oscil¬ 
lators. Note that the real part of the complex conjugate 
eigenvalues /ii ^2 is positive and smaller than 1 and in fact 
slowly goes to 0 as a approaches the stabilizing Hopf bi¬ 
furcation, see Fig. where Re jxi is depicted in green 
(dashed) and crosses the real axis at the Hopf bifurca¬ 
tion point. 

The (slow) unstable eigenvalue /ii is compensated by 
the (fast) stable eigenvalue /is. The transient time T ^ 
1/Re/ii is shorter if the real part of fii increases. If 
Re /ii goes to zero (which it does at the Hopf bifurcation 
point), the transient time goes to infinity. In contrast, 
the transient time increases with the absolute value of 
the stable eigenvalue. 

We considered an additive Gaussian randomization of 
initial data, with fixed variance. Numerical experiments 
then indicated that the observed mean transient time is 
proportional to the ratio of the real parts of the eigenval¬ 
ues, more precisely, the larger the ratio, the longer the 
transient time. The stochastic dependence of T on the 
initial conditions, i.e., the position and variance of the 
distribution of the oscillators will not be pursued in this 
publication. We include this by a factor /d > 0, which 
so far we have determined numerically, by fitting one or 
more values of a. Here f3 depends only on the initial 
conditions and is independent of P, n, N and a. 

We conclude that the dependence on P and a of the 
expection value (T) of the transient time T can be de¬ 


scribed by 


(T) = /3 IRe/isl/Re/Ui 


= /3 


2 



— (jO 


2 


—A T 2 k,p yi(t — 


2 



( 20 ) 

( 21 ) 


Note that, due to the sensitive dependance on the ini¬ 
tial conditions, we only give an estimate for the expec¬ 
tation value of the transient time T. Fig. shows the 
numerically calculated transient times vs a for different 
variances of the initial distribution. The solid curves rep¬ 
resent Eq. (21) with fit parameters (3. 


VII. CONCLUSION 

We have investigated a network of nonlocally coupled 
Stuart-Landau oscillators under a coupling which breaks 
-symmetry. Following numerical simulations with clus¬ 
tered initial conditions, we have found a family of oscil¬ 
lation death states (inhomogeneous steady states) which 
can be distinguished by their cluster size. Clustering al¬ 
ways emerges in pairs of in-phase and anti-phase clusters. 
Two different cases occur: Either the pattern of the ini¬ 
tial cluster configuration coincides with the asymptotic 
cluster distributions, or additional clusters appear in the 
middle of the original clusters. The shape of the clus¬ 
ters depends on the coupling range P, i.e., the number 
of nearest neighbors which are coupled. If P is either an 












integer multiple of the cluster size, or an integer multiple 
minus one, then all oscillators within the cluster exhibit 
the same radius and phase. In all other cases a linear 
cluster deformation (branch splitting) occurs. We have 
developed an approximate analytical description by a re¬ 
duced model of two mean-field coupled Stuart-Landau 
oscillators. In particular, we have analytically calculated 
the onset of stable oscillation death with prescribed clus¬ 
tering. Our theory goes beyond standard mean-field the¬ 
ory, which only gives a rough estimate. Specifically, we 
have extended the mean-field coupling by an approxima¬ 
tion for the spatial deviation from the mean. As a con¬ 
sequence, we are able to predict the boundaries of the 
different stability regimes of m-cluster oscillation death 
analytically with high precision. 

In addition to the asymptotically stable oscillation 
death states, we also find a region where oscillation death 
is transient: It persists for a long time but then disap¬ 
pears in favor of synchronized oscillations. The tran¬ 
sient behavior occurs due to the interaction of a slow 
unstable eigenvalue and a fast stable eigenvalue of the 
inhomogeneous steady state. We have analytically cal¬ 
culated a scenario leading from the transient behavior to 
the asymptotically stable oscillation death with the cou¬ 
pling strength serving as a bifurcation parameter. Near 
the stabilizing Hopf bifurcation, the calculated transient 
time goes to infinity which is in accordance with the nu¬ 
merical simulations. 

Note that in this paper we have confined attention 
to spatially coherent clusters of inhomogeneous steady 


states. As shown in [nmi], there also exist hybrid states 
consisting of coexisting domains of spatially coherent and 
incoherent oscillation death, called chimera deaths for cer¬ 
tain initial conditions. This is an indication of the high 
multistability of the system. The regimes of existence of 
1-cluster, 3-cluster, 5-cluster etc. chimera death in the 
(cr, P) parameter plane is similar to Fig. In addition, 
for small P there exist transient amplitude chimeras of 
coexisting domains of spatially coherent (synchronized) 
and incoherent oscillations. 

Stuart-Landau oscillators as normal form of systems 
near Hopf bifurcation can serve as a model system for 
many applications. The set-up for both the mean-field 
theory and its correction does not depend on the spe¬ 
cific form of Stuart-Landau oscillators. Similar symmet¬ 
ric cluster states, maybe even including more branches, 
can be expected in any network with a clearly defined 
symmetry, in which case the methods presented in this 
paper can be employed. 
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